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Abstract. We present ISO mid-infrared spectrophotometry and far-infrared photometry of galaxies with Hidden 
Broad Line Regions (HBLR). We also present radiative transfer models of their spectral energy distributions 
which enable us to separate the contributions from the dusty disc of the AGN and the dusty starbursts. We find 
that the combination of tapered discs (discs whose thickness increases with distance from the central source in the 
inner part but stays constant in the outer part) and starbursts provide good fits to the data. The tapered discs 
dominate in the mid-infrared part of the spectrum and the starbursts in the far-infrared. After correcting the 
AGN luminosity for anisotropic emission we find that the ratio of the AGN luminosity to the starburst luminosity, 
LAGN/L sb , ranges from about unity for IRAS 14454-4343 to about 13 for IRAS01475-0740. Our results suggest 
that the warm IRAS colours of HBLR are due to the relatively high Lagn / L s b- Our fits are consistent with the 
unified model and the idea that the infrared emission of AGN is dominated by a dusty disc in the mid-infrared 
and starbursts in the far-infrared. 
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1. Introduction 

Since the discovery of broad lines in the polarized flux of 
the prototypical Seyfert 2 galaxy NGC1068 by Antonucci 
, & Miller (|1985[) much effort has been directed towards 
the study of this phenomenon in large samples of narrow- 
lined A GN (Miller & Goodrich ITMH Young et al. llj)i)fJal 
I1996bl Heisler, Lumsden & Bailev HWH Tran 127711111217051 
Lumsden et al 200TJ . The main aim was to test the unified 
model of active galaxies according to which many of the 
differences between narrow-lined and broad-lined AGN 
can be understood as arising from the presence of an ob- 
scuring torus viewed from different directions (Antonucci 

The presence of a geometrically and optically thick 
dusty torus implies that a large fraction of the optical and 
ultraviolet radiation emitted by the central engine of the 
AGN is absorbed by dust and re-radiated anisotropically 
in the infrared (Efstathiou & Rowan- Robinson 1990, 1995 
Pier & Krolik 1TM21 Granato & Danese \TM3$. Detailed 
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radiative transfer calculations show that the predicted in- 
frared spectrum of a dusty torus is sensitive to a number 
of parameters including the inclination. Deep absorption 
features at around 10 \im due to silicate dust are predicted 
for edge-on views of the torus. For face-on views a variety 
of spectra are predicted ranging from emission/absorption 
features at lO/im to completely featureless spectra. 

The spectral energy distributions (SEDs) predicted by 
most radiative transfer models of dusty discs and tori pre- 
sented to date are generally narrower than the observed 
infrared SEDs. This is not necessarily a limitation of the 
models. NGC1068 has in fact been known for some time 
to be a composite object (e.g. Telesco et al. I1984fl where 
the Seyfert nucleus is surrounded by a 3 kpc diameter 
ring of star formation. The AGN dominates in the mid- 
infrared part of the spectrum (Rieke & Low 1975) and the 
starburst dominates in the far-infrared. 

Rowan- Robinson & Crawford (1989) showed that a 
three component model (cirrus, starburst and AGN) could 
explain observations of IRAS galaxies. Genzel et al. ([1998) 
argued on the basis of ISO spectroscopy that about half 
of the ultraluminous infrared galaxies harbour simultane- 
ously an AGN and starburst activity. ISO spectroscopy 
(see also Roche et al. I1991fl also showed the prevalence 
of mid-infrared emission features due to PAHs in AGN 
(Clavel et al.l2TT00). These features are believed to be emit- 
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ted by the starburst in these systems as the radiation field 
is too hard for the PAH molecules to survive in the dusty 
torus (Roche et al. 119911 Siebenmorgen et al. 12004ajl . 

The coexistence of dusty AGN and nuclear starbursts 
and the low spatial resolution of current infrared obser- 
vations require a composite star burst/ AGN torus model 
to explain the observed SEDs. In recent years there has 
been considerable progress in attempts to model the 
SEDs emitted by starburst galaxies (Rowan-Robinson & 
Efstathiou 119981 Kriigel & Siebenmorgen 1994; Silva et 
al. 119981 Efstathiou, Rowan-Robinson & Siebenmorgen 
2000 Siebenmorgen, Kriigel & Laureijs 1200 lfl and axi- 
symmetric dust distributions in AGN (Pier & Krolik 1992 
Granato & Danese 119941 Efstathiou & Rowan- Robinson 
1995J. These models have been successfully compared to 
the SEDs of a number of infrared galaxies. 

The results from the different polarimetric surveys 
show that not all narrow-lined AGN show broad lines 
in polarized flux. Galaxies known to have HBLR gener- 
ally show warmer infrared colours compared to galaxies 
without HBLR. The first interpretation of this result was 
that the occurrence of the phenomenon requires a special 
viewing geometry (Heisler et al. 11997}. Later Alexander 
(HU1), Lumsden et al. lj2TjT]T)l and Tran lj2lMl |2TM|> ar- 
gued against this model. The cold colours of galaxies with- 
out HBLR were instead mainly attributed to the domi- 
nance of the host galaxy. Galaxies without HBLR have in 
general weaker AGN. This is consistent with the compos- 
ite starburst/ AGN torus model. Tran went a step further 
and argued that narrow-lined objects without HBLR do 
not have a broad line region at all and therefore the unified 
model does not apply to all AGN. 

In this paper we report Infrared Space Observatory 
(ISO) mid-infrared spectrophotometric and far-IR pho- 
tometric observations of two galaxies (IRAS 14454-4343 
and 3C321) that show HBLR. We use radiative transfer 
models to separate the contributions from starbursts and 
AGN to the SEDs of these two galaxies. We also model 
IRAS01475-0740 for which we extracted spectrophotomet- 
ric data from the ISO archive. The data presented here 
allow us to carry out the most stringent test to date to 
the composite starburst /AGN torus model for the infrared 
emission of AGN. We assume a flat Universe with A = 0.73 
and H = 71Km/s/Mpc. 

2. Observations 

Mid infrared spectrophotometric imaging of 
IRAS 14454-4343 and 3C321 were performed on 03- 
Sep-1997 using the ISOCAM circular variable filter 
(CVF, Cesarsky et al. HMty . The ISO archive number 
(TDT) of the observations are 65800208 for 3C321 
and 65800106 for IRAS14454-4343, respectively. The 
observations cover a total field of 96" x 96". For each CVF 
step between 16.3/mi down to 6.4/im about 30 exposures 
of 2.1s integration time using the 3" lens were read out. 

The ISOCAM data were reduced with the ISOCAM 
Interactive Analysis (CIA, Ott et al. I1997[> . We used the 



default data reduction steps of CIA: dark current subtrac- 
tion, initial removal of cosmic ray hits (glitches), detector 
transient fitting, exposure coaddition and flat fielding. 

The dark current depends on the orbital position of 
the ISO space-craft and the temperature of the ISOCAM 
detector. The applied correction is based on a model de- 
scribed by Roman & Ott (1999). The deglitching is done 
by following the temporal signal variation of a pixel using 
a multi-resolution wavelet transform algorithm (Starck et 
al. H997fl . The response of the detector pixels strongly de- 
pends on the previous observations and there are long 
term hysteresis effects for each detector clement after 
changes of the photon flux level. The detector flux tran- 
sient fitting method for ISOCAM data was developed by 
Coulais & Abergel (2000). After application of the de- 
fault deglitcher some residuals of cosmic ray impacts were 
still visible in the data. Therefore after the detector flux 
transient correction we applied a second cosmic ray rejec- 
tion method which is basically a multi-sigma clipping of 
the temporal signal (Ott et al. I2000f) . To determine the 
source spectrum, we perform a multi-aperture photom- 
etry on each image of the different CVF steps. For all 
apertures centered on the brightest pixel of the source we 
determine the background as the mean flux derived in a 
4 pixel wide annulus which is put 2 pixels away from the 
greatest aperture. In this way, we obtain aperture fluxes 
that flatten with increasing aperture radius and approach 
an asymptotic value. The same procedure is repeated on a 
theoretical and normalised point source image calculated 
for the central wavelength of each CVF image (Okumura 
1998). This provides a correction factor of beam effects of 
the multi-aperture analysis. In order to derive the statis- 
tical uncertainty of the flux, we perform the same proce- 
dure on the rms image and quadratically add the rms of 
the background estimate. The so far reduced flux must be 
converted into astronomical units (mJy) by multiplication 
with the standard conversion factor (as given in the cali- 
bration files: CCGLWCVF1 and CCGLWCVF2 provided 
within CIA; see ISOCAM Handbook, Blommaert et al. 

(Ennui. 

Both targets were also observed with the PHT in- 
strument at 170/iTO in mini-raster mode. The TDT of 
the observations was 61801012 for IRAS14454-4343 and 
60201514 for 3C321. IRAS14454-4343 was also observed at 
65 fim (TDT 61801013) and 3C321 at 90/xm (60201515). 

For the ISOPHT data standard PIA version 9 data re- 
duction was used. The "actual" FCS responsivity, and not 
the default average one, was used. The photometry was 
derived from the small maps assuming unresolved point 
sources. There was no indication for a significant extended 
resolved flux contribution in the maps. The results of the 
ISOPHT data reduction are given in Table 1. 



A. Efstathiou and R. Siebenmorgen: ISO Observations and models of HBLR 



3 



Table 1. ISOPHT photometry of IRAS 14454-4343 and 
3C321. 



Name 




<S90/Jm/Jy 


Si70 Mm /Jy 


IRAS 14454-4343 


4.5 ±0.3 




2.3 ±0.35 


3C321 




0.78 ±0.23 


0.4 ±0.12 



3. Radiative transfer models 

3.1. Starburst models 

Efstathiou, Rowan- Robinson & Siebenmorgen ( 2000 here- 
after ERS00) developed a starburst model that includes 
stellar population synthesis (Bruzual & Chariot Il993|) . a 
simple model for the evolution of giant molecular clouds 
and detailed radiative transfer that includes the effect of 
transiently heated particles/PAHs. The models make pre- 
dictions about the emitted spectrum of starburst galaxies 
from the UV to the submillimeter as a function of the star 
formation history of the starburst and its age. Some illus- 
trative results are given in Fig. 1. The effect of age on the 
SED is to make it progressively cooler with the interest- 
ing prediction that the submillimeter luminosity becomes 
almost independent of starburst age. Note also that the 
absorption feature at 10/ito gets shallower with age. The 
models assume that the dust to gas ratio does not evolve 
with age. 

The ERS00 models were applied with considerable suc- 
cess to fit the SEDs of the nearby starbursts M82 and 
NGC6090. The models were also used to model the SEDs 
of ultraluminous (Farrah et al. 12003)1 and hyperluminous 
(Rowan-Robinson I2UUU1 Verma et al. I2TJU21 Farrah et al. 
2002) infrared galaxies. In this study we consider models 
of the type described in ERS00 with a range of average 
molecular cloud densities (n av — 300 — 1500cto~ 3 ) and 
molecular cloud masses (M GM c = 10 7 - 10 8 M Q ). The 
star formation efficiency r/ is kept fixed at 0.25. The mod- 
els assume two more parameters: the time constant of the 
exponentially decaying star formation rate (t = 1 — 30 
Myrs) and the age of the starburst (t = — 72Myrs). 

3.2. AGN models 

Models of the infrared emission from dusty tori or discs 
in active galactic nuclei were presented by Pier & Krolik 
iflU^l lllMJl . Granato & Danese IT^Ml Efstathiou & 
Rowan- Robinson ll9951 van Bemmel & Dullcmond I2003f) . 

One of the principal constraints for these models 
were spectrophotometric observations of AGN at 8-13^m 
(Roche et al. H991f) . These observations showed moderate 
absorption features at 10/im in type 2 AGN, which are at- 
tributed to silicates, but rather featureless spectra in type 
1 AGN. 



Pier & Krolik (1992) presented models of moderate 
to high optical depth tori which were approximated as 
cylinders of uniform density. The most optically thick of 
their models produced flat 8-13/im spectra or moderate 
absorption features when the tori were viewed face-on but 
the SEDs were rather narrow. In their application of these 
models to NGC1068, Pier & Krolik ip33|) introduced an 
additional optically thin component that broadened the 
emission from the optically thick torus. 

Granato & Danese 1)1994)1 presented models for flared 
discs of moderate optical depth. With a standard inter- 
stellar dust mixture such models predict emission features 
at 10/ito when the discs are viewed face-on (Efstathiou 
1990). Granato & Danese showed that these features are 
suppressed if the silicates are destroyed (for example by 
shocks) in the inner part of the dusty disc. Grain mixtures 
that differ from the standard galactic one have also been 
discussed by Laor & Draine 1)19931) . Maiolino et al. H2001al 
I2001bj) . van Bemmel & Dullemond l|2TiU3"jl and Galliano et 
al. 1)2003)1 . Granato, Danese & Franceschini (1994) pre- 
sented a fit to the SED of the nucleus of NGC1068 with a 
model that assumes that the dust density varies exponen- 
tially with a function of the polar angle O. 

Efstathiou & Rowan-Robinson 1)19951 hereafter 
ERR95) considered anisotropic spheres (see also Rowan- 
Robinson et al. i995J, flared and tapered discs. Flared 
discs have a thickness that increases linearly with distance 
from the central source. The thickness of tapered discs 
also increases with distance in the inner disc but assumes 
a constant value in the outer disc. ERR95 concluded that 
the geometry that best fits the observational constraints 
is that of tapered discs. Efstathiou, Hough & Young 
( 1995) proposed a model for the nucleus of NGC1068 that 
involved a tapered disc and dust in the ionization cones. 
The presence of conical dust in NGC1068 is supported by 
subarcsecond mid-infrared imaging (e.g. Cameron et al 
1993). If the conical dust is diffuse and composed of the 
standard interstellar mixture it will emit a strong silicate 
emission feature. This dust should therefore either be 
concentrated in clouds or be composed of large grains 
only. Efstathiou, Hough & Young attempted to model 
approximately the emission of conical dust by increasing 
the abundance of very large (30/ito) grains in the cones. 

Rowan-Robinson (199SJ) presented a model for the in- 
frared emission of PG quasars that treated an anisotropic 
distribution of clouds which were approximated by nar- 
row shells. This model also produces weak silicate emis- 
sion features and broad SEDs. The idea of a clumpy torus 
was explored in much more detail by Nenkova, Ivezic & 
Elitzur (2002} who calculated the emission of individual 
clouds with the slab approximation. The latter authors 
also computed the emission/absorption by an ensemble of 
clouds in a flared disc. 

The AGN models considered in this study are the ta- 
pered disc models described in ERR95 that represent one 
possible solution to the problem. For simplicity we do not 
consider conical dust in this paper. The reader is referred 
to ERR95 for a discussion of how the predicted SEDs de- 
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Fig. 1. Spectral energy distributions of Efstathiou et al. starbursts in which the star formation rate decays exponen- 
tially with time (time constant 30Myrs). The SED is shown for three different starburst ages: lOMyrs (solid), 40Myrs 
(dotted) and 70Myrs (dashed). The models assume a molecular cloud mass Mgmc of 2 x 10 7 M©, a star formation 
efficiency rj of 0.25 and an average molecular cloud density fi av of 600cm 3 . The initial star formation rate is assumed 
to be 100 M Q /yr. 



pend on the assumed parameters. We consider tapered 
disc models in which we vary two parameters: the half- 
opening angle O which we vary in the range 30°, to 60° 
(this is related to the angle Oi defined by ERR95 by 
6 = 90° — 9i) and the inclination i which covers the 
range to tt/2 (face- to edge-on tori). The SED is com- 
puted for 37 lines of sight equally spaced in i. The ratio of 
inner to outer disc radii T\jri is fixed at 0.05 which is the 
middle of the values considered by ERR95. The equatorial 
1000A optical depth is kept fixed at 1000 and the ratio 
of disc height to inner radius h/r\ is assumed to be 10 for 
all models. As in ERR95 the density distribution in the 
tapered disc is assumed to follow r . 



Galliano et al. f2003) carried out a comprehensive 
study of the prototypical Seyfert 2 galaxy and HBLR 
NGC1068. They found that the nuclear SED can be fit- 
ted successfully with both tapered discs and flared discs 
whose density varies with polar angle 0. Galliano et al. 
found that the size and optical thickness of the torus are 
reasonably well constrained by the data but the geome- 
try of the torus is not. Good fits are obtained with models 
that assume that the inclination is in the range 30° to 50°. 
For reference, the model of Efstathiou, Hough & Young 
(1995J) for the nucleus of NGC1068 assumed an inclina- 
tion of 45°. Because we are only considering a small set of 
tapered disc models in this paper, our predicted inclina- 
tions may be uncertain by 10-20°. This also introduces an 
uncertainty in the anisotropy corrected luminosity of the 
AGN (see section 4.1). 



4. Results of model fitting 

Given the set of starburst and torus models described 
above we seek the combination which provides the min- 
imum reduced \ 2 - The torus emission can vary with 
the inclination but the starburst emission is assumed to 
be inclination- independent. Furthermore, we assume that 
there is no absorption of radiation emitted by the torus 
by the starburst and vice versa. This is not an important 
limitation if the starburst is concentrated in a circumnu- 
clear ring as in the case of NGC1068, although in such 
a case we may expect some heating of the outer parts of 
the torus by the starburst. If, however, a significant frac- 
tion of the radiation emitted by the AGN in directions 
not covered by the torus is intercepted by circumnuclear 
or galactic dust then our estimated starburst luminosity 
will be overestimated. 

To be consistent with the fact that no broad lines are 
seen directly in these objects, we only consider models 
where the line of sight to the central source passes through 
the disc. The resulting fits are given in Fig. 2 and the best 
fit parameters are tabulated in Table 2. There are no other 
combinations of starburst and AGN torus/inclination that 
give a x 2 lower than Xmin + 1- The best fit models are 
found to be in good agreement with the broad band SEDs 
as well as the detailed ISO spectrophotometry in the 5- 
13/im range. 

4.1. IRAS 14454-4343 

Polarized Broad H a and Hp were discovered in 
IRAS14454-4343 by Young et al. l|1996ajl . A very good fit 
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log (l/iun) We" 

Fig. 2. Best fits to the rest frame SEDs of IRAS 14454-4343, 3C321 and IRAS01475-0740 with a starburst/AGN 
torus combination model. The total emission is plotted with a solid line whereas the individual contributions from the 
starburst and AGN are given by the dashed and dotted lines respectively. In the case of 3C321 and IRAS01475-0740 
where optical/near-IR data is available, the emission by an old stellar population (12.5Gyrs) normalized to the H band 
is shown with a dash-dot line. 
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Table 2. Best fitted AGN torus parameters and derived luminosities: A is the anisotropy factor (see section 4.1), i 
is the angle between the line of sight and the polar axis of the disc, Q is the half-opening angle, L s {,, Lagn and L to t 
are the starburst, AGN and total 1-1000/im luminosities respectively. Note that the emission from the AGN torus 
is anisotropy corrected (see section 4). For the starburst models the best-fitted parameters for all objects are: time 
constant of the exponentially decaying star formation rate is 3 x 10 7 years, age is 7.2 x 10 7 , average molecular cloud 
density 600cm~ 3 and molecular cloud mass 2 x 10 7 M Q . 



Name 


z 


Xmin /^-f 


.4 


i 


e 


log Lagn/L & 


log L 3b /L e 


log Ltot/L® 


IRAS 14454-4343 


0.03856 


1.19 


1.4 


58° 


45° 


11.40 


11.37 


11.69 


3C321 


0.09610 


1.26 


7.3 


73° 


60° 


12.30 


11.40 


12.35 


IRAS01475-0740 


0.01766 


0.85 


1.7 


60° 


45° 


10.78 


9.67 


10.81 



to the spectrophotometry and far-IR photometry is ob- 
tained with a model where i is 58°. 

As it is clear from the plot of the spectral energy dis- 
tribution of the torus for different inclinations given by 
ERR95, the emission from the torus is highly anisotropic. 
We quantify this with the anisotropy factor A(i) which is 
defined to be 



2 fo /2 S(i) smi di 
irS(i) 



A(i) 



where S(i) is the bolometric emission at angle i. 

For the specific model and inclination that gives the 
best fit to IRAS1445-4343, A is 1.4 which implies that the 
AGN is about 40% more luminous than it would be if the 
observed emission was radiated isotropically. After this 
correction the AGN torus is found to be approximately 
as luminous as the starburst. The predicted intrinsic lu- 
minosity of IRAS14454-4343 (4.9 x 10 n L Q ) places this 
object in the Luminous Infrared Galaxy (LIRG) class. 

The starburst age is predicted to be 72Myrs and the 
initial visual optical depth of the molecular clouds is pre- 
dicted to be 250. The latter is a factor of five higher than 
the optical depth of molecular clouds in the M82 and 
NGC6090 starburst (ERS00). More details of assumed pa- 
rameters and derived luminosities are given in Table 2. 

4.2. 3C321 

A polarized broad component of H a was discovered in 
this powerful FR II galaxy by Young et al. I|1996bjl . The 
spectrophotometric data for this object are not as good 
as for IRAS1445-4343 as this object is about an order of 
magnitude fainter in the mid-infrared. Nevertheless a clear 
10/zm feature in absorption is detected. The best fit model 
is the most geometrically thin of the models considered i.e. 
the model with o = 60°. The best fit i in this object is 
73°. This implies a higher anisotropy factor which is found 
to be 7.3 in this case. After correcting for this effect we find 
that the AGN is about a factor of 8 more luminous than 
the starburst. The predicted intrinsic luminosity of 3C321 



(2.2 x 1O 12 L ) places this object in the Ultraluminous 
Infrared Galaxy (ULIRG) class. 

The best fit starburst parameters for 3C321 are iden- 
tical to those found for IRAS14454-4343. 

4.3. IRAS01475-0740 

Polarized broad lines in this object were discovered by 
Tran ( 200.1). The best fit disc model is identical to that 
for IRAS1445-4343 but i is slightly higher (60°) i.e. this 
object is viewed more edge-on. The anisotropy correction 
factor is also slightly higher (1.7). This object is the least 
luminous of the objects studied with a total intrinsic lumi- 
nosity of 6.5 x 10 10 Lq. Interestingly, however, it has the 
highest ratio of AGN to starburst luminosity (13). The 
fitted starburst parameters are also identical to those for 
IRAS14454-4343 and 3C321. 

5. Discussion and conclusions 

We presented ISO mid-infrared spectrophotometry and 
far-IR photometry for two objects (IRAS14454-4343 and 
3C321) that show (polarized) Hidden Broad Line Regions. 
We also extracted from the ISO archive mid-IR spec- 
trophotometry for a third object (IRAS01475-0740) that 
shows the same phenomenon. We used radiative trans- 
fer models of torus and starburst emission to separate 
their contribution to total luminosity. We find that the 
combination of the tapered discs of Efstathiou & Rowan- 
Robinson l|1995[) and the starburst models of Efstathiou, 
Rowan-Robinson & Siebenmorgen (2000) provide good 
fits to the data. These results support the idea of a 
two-component model to explain the infrared emission of 
AGN. The same model has been applied to fit the SEDs 
of ultraluminous (Farrah et al. I2003f) and hyperluminous 
(Morel et al. I2TMTT1 Farrah et al 12711721 Verma et al HEJJ) 
infrared galaxies although in those studies mid-infrared 
spectrophotometry was not available. A similar study was 
carried out by Rowan-Robinson (1995) who modeled the 
SEDs of PG quasars with a combination of the starburst 
model of Rowan-Robinson & Efstathiou ( 1993J) and a 
model for the AGN emission that considered an ensemble 
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of geometricaly thin shells. Models of the type presented 
here were also applied to fit the SEDs of Cen A (Alexander 
et al.EiSJl and Circinus (Ruiz et al. l5tM|l . 

Nenkova et al. H2002JI presented a model that is based 
on the idea that the dust in AGN is concentrated in clouds. 
They found that if the dust is contained in 5-10 clouds 
(each with a visual optical depth greater than 60) along 
radial rays through the torus they can produce SEDs ex- 
tending to the far-infrared. The broadest of the predicted 
SEDs shows a silicate feature in emission for type I ob- 
jects. More recently, Dullemond & van Bemmel (2005) 
presented models of clumpy tori around active galactic 
nuclei and compared them with equivalent smooth tori. 
They found that the lOfim emission feature of the clumpy 
models, when seen almost face-on, is not appreciably re- 
duced compared to the equivalent smooth models. 

Kuraszkiewicz et al. (2003) presented flared disc mod- 
els for three X-ray selected active galaxies. The mod- 
els reproduce the whole infrared SED without the need 
for a starburst component. The models presented by 
Kuraszkiewicz et al. show strong silicate features in emis- 
sion. The latter authors note that these features would 
be eliminated if the torus were clumpy but as discussed 
above this may not necessarily be the case. 

Siebenmorgen et al. I|2004bp presented radiative trans- 
fer models for 3CR radio galaxies detected by ISOCAM. 
Although the models of Siebenmorgen et al. are spheri- 
cally symmetric, they take into account for the first time 
the destruction of small grains and PAHs by the radia- 
tion field of the central engine of the AGN. Their models 
show that when the dust is at large distances from the 
central engine (several kpc) the SED peaks in the far-IR 
and strong PAH features are predicted. When the dust is 
close to the central engine, a situation that approximates 
a nuclear torus, the SEDs shift to shorter wavelengths and 
do not show strong PAH bands. 

It would be interesting to see if the models considered 
in this paper and those discussed above can successfully 
model the detailed SEDs of AGN that will soon become 
available with the Spitzer Space Telescope. 

The results from the various spectropolarimetric sur- 
veys carried out to date show that not all narrow-lined 
AGN contain HBLR. Heisler et al. i)l»97|l pointed out that 
galaxies that harbour HBLR tend to have warm mid- to 
far- infrared colours (i^o /-F25 < 4) . To interpret this result 
they proposed a model according to which the orientation 
of the dusty torus is the main factor that determines the 
visibility of polarized broad lines. Heisler et al. suggested 
that the reason for this is that scattering takes place very 
close to the nucleus. In their model non-HBLRs are those 
narrow-lined AGN which are viewed at high inclination 
and therefore the scattered (and hence polarized) radia- 
tion is also blocked by the torus. In a subsequent paper 
Lumsden et al. H2001[) presented results from a spectropo- 
larimetric survey of a complete far-infrared selected sam- 
ple of Seyfert 2 galaxies. This study confirmed the findings 
of Heisler et al. that galaxies with HBLR tend to have 
warm mid- to far-infrared colours. However, Lumsden et 



al. presented evidence showing that the warmth of the in- 
frared emission is not solely due to inclination. They found 
instead that the main determining factor is the relative lu- 
minosity of the AGN and the host galaxy as proposed first 
by Alexander l|2001[) . HBLRs are found in objects where 
the AGN dominates the infrared luminosity. Lumsden et 
al. also found that HBLRs have less obscuration than non- 
HBLRs of the same luminosity suggesting, that orienta- 
tion may also be important if obscuration arises in a torus. 

Recently Lumsden, Alexander & Hough H2004|) pre- 
sented results from a spectropolarimetric survey of a sam- 
ple of Compton-thin Seyfert 2 galaxies. They found that 
the detection of HBLR in this sample is considerably 
higher than in Seyfert 2 galaxies as a whole. They also 
found HBLR in galaxies with cool mid- to far-infrared 
colours {F eo /F25 > 4). These observations argue against 
the idea proposed by Tran ( 2003 ) that there are two pop- 
ulations of Seyfert 2s, those with a broad line region and 
those without one. 

The available evidence therefore suggests that the oc- 
currence of the HBLR phenomenon requires a strong AGN 
relative to the host galaxy and possibly lower obscura- 
tion which may imply a dependence on inclination. The 
method applied in this paper allows a determination of 
both the ratio of AGN to starburst luminosities and the 
inclination. So a similar analysis on a larger sample of 
Seyfert 2s preferably with data provided by new space 
telescopes, such as the Spitzer Space Telescope, promises 
to be useful for understanding the HBLR phenomenon. 

The tendency of galaxies with HBLR to show both 
warm mid- to far-infrared colours and evidence for strong 
AGN relative to the host galaxy (as measured by the 
[OIII] 5007 A to H^ ratio) supports the idea of the 
two-component model for the infrared emission of AGN. 
Galaxies with warm colours are objects with a strong AGN 
whose emission peaks in the mid-infrared and therefore 
boosts the 25/ito flux. 

As it is evident from Fig. 2, even in objects in which the 
AGN dominates the bolometric luminosity, the starburst 
dominates at wavelengths longer than 60-100/im. By con- 
trast the rest-frame mid-infrared emission is dominated 
by the AGN torus. If we assume that the same model ap- 
plies to the high redshift Universe, this result justifies the 
use of submillimctcr and millimeter fluxes to derive star- 
formation rates for high redshift objects found in blank 
field surveys (e.g. Hughes et al. lEMl Scott et al.EIEI) or 
pointed observations of high redshift objects (e.g. Priddey 
& McMahonEEH)- 
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